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We apply a renormalization group approach to the determination of the phase diagram of the 
t — t' Hubbard model at the Van Hove filling, as function of t' /t, for small values of U/t. The model 
presents ferromagnetic, antiferromagnetic and d-wave superconducting phases. Antiferromagnetism 
and d-wave superconductivity arise from the same interactions, and compete in the same region of 
parameter space. 
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In receat years great effort has been devoted to study the role of Van Hove singularities in two-dimensional electron 
liquidstrEl. Most part of the interest stems from the evidence, gathered from photoemission exp|eriments, that the 



I . hole-doped copper oxide superconductors tend to develop very flat bands near the Fermi levelQu. The model most 
widely used for the Cu-0 planes is the Hubbard model, which also shows strong antiferromagnetic correlations. The 
^ , correspondence with the phenomenology of the high-Tc materials is supposed to be attained in the intermediate 
' to strong coupling regime. Weak coupling approaches to the model hai|e shown that it is more likely to develop a 

a] spin-density-wave instability than superconductivity near half-fillin^'ElEil. 
The inspection of the actual Fermi line of most of the copper oxide superconductors shows, however, a sensible 
departure from nesting. In the context of models with on-site interaction, the dispersion relation seems to be best 
^ fitted by introcLucing both nearest neighbor t and next-to-nearest neighbor t' , accomodated in the so-called t — t' 
O Hubbard modellj. In this model, the absence of perfect nesting implies that the antiferromagnetic correlations have 
^ O^ , a less divergent behavior when the Fermi level lies at the Van Hove singularity. For this reason, it is best suited for 
the implementation of a renormalization group (RG) approach to the ground state ppaperties, that is needed anyhow 
T-H to deal with the logarithmic singularities arising from the divergent density of stateal3. 

J> ■ In this Letter we look for the instabilities of the t — t' Hubbard model, filled up to the level of the Van Hove 
in singularity, following the wilsonian RG program of Refs. [l^ , |l5| . Our approach is that of integrating virtual states of 
. two energy slices above and below the Fermi level in an energy range given by the cutoff Ec, Ec ~ \dEc\ < \E\ < Ec. 
' We are interested in the scaling behavior of the interactions and correlations under a progressive reduction of the 
. cutoff, which leads to the description of the low-energy physics about the Fermi level. 

' In the RG process, we have to make sure first that the interaction does not scale to zero at the classical level, that 
^\ , is, it displays marginal behavior as we approach the Fermi line. When the cutoff is sufficiently close to the Fermi 
■ level, as shown in Fig. ^ most part of the states at hand are in the neighborhood of any of the two Van Hove points 
(^lO) (^) and (0,7r) (B), and the dispersion relation may be approximated by two respective patches about each of 
them: 

^\ eA.B{k)^TitT2t')kla^±it±2t')ky (1) 



where a is the lattice constant. The angle between the two separatrices is 2 arctan[(t + 2t') /{t — 21')]. The low-energy 
effective action may be written in the form 



/ dcjd^fc^ (w a+^(k,a;)aa,a(k,a;) - £a(k) a+^(k,u;)aa,<T(k,a;)) 

-U I duj(fkp'i{k,uj)pi{-k,-uj) (2) 



where aa^a-{o-a a) ^.n electron annihilation (creation) operator {a labels the Van Hove point), p-\^i are the density 
operators in momentum space, and the momentum integrals are restricted to modes within the energy cutoff, \ea (k) | < 
Ec. It is a fact that the action (J^ defines a model that is scale invariant at the glassical level, that is, the different 
terms that build up (||) are invariant altogether under a rescaling of the cutoff E^. This implies that the interaction 
in the t — t' Hubbard model has marginal behavior, at least at the classical level, and it is therefore susceptible of a 
low-energy description within the RG framework. 
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In our low-energy model, the two patches for the dispersion relation give rise to an additional flavor index. The 
number of possible interactions proliferates accordingly. The different terms are depicted in Fig. ^ . It might seem 
that the analysis of the quantum corrections in the model should bear a great similarity with the g-ology description 
of one-dimensional models. In two dimensions, however, kinematical constraints play a crucial role. Performing the 
wilsonian RG analysis, no singular cutoff dependences are found in the particle-particle channel, for instance, unless 
the two colliding particles have opposite momenta. This conclusion may be drawn from the general discussion carried 
out in Ref. The BCS instability is to be analyzed, then, as a singularity in a response function computed at a 
definite value of the momentum. On the other hand, singular dependences on the cutoff are found in the particle-hole 
channel, when the momentum transfer is about zero or about the momentum connecting the two Van Hove points. 

The interactions shown in Fig. ^ are renormalized by the high-energy modes at the quantum level. In all the 
four cases, the only contribution to second order in perturbation theory which is ~ dEc/Ec is given by a diagram of 
the type shown in Fig. ^. As mentioned before, the contributions in the particle-particle channel are 0[[dEcf') for 
generic values of the incoming and outgoing momenta, and there arc no more diagrams which can be built up from 
the interaction in (||). To second order we get the RG flow equations 



d Uintras _ 1 fjj^2 , jT^ 

dEc ~ 2'kH" 

dUback ^ ^TT TT \ ^ A\ 

OEc TT^t 

„ dUinters _ I , /,.2 , rr2 \ 
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are the prefactors of the polarizabilities at zero 



where c = l/^l -A{t'/tY and c' = log + ^1 - A{t' /tf^ li^t'/t) 

and Q momentum transfer, respectively. 

The RG equations (^)-(|^) describe a flow that drives the couplings to large values, as the cutoff is sent to the 
Fermi line. This is due to the localization in space of the interactions, which makes them strongly spin-dependent. 
The RG flow of extended, spin- independent, interactions is towards lower couplings, as discussed in Ref. 0. In that 
case, RPA-like screening is the dominant effect. The limit of a short range interaction is special in that there is 
a cancellation between the 'particlc-hole bubble' and 'vertex-correction' second order diagrams. For a purely local 
interaction like that of the t — t' Hubbard model, the only contribution left correspons to the antiscreening diagram 
shown in Fig. ^. 

In order to apply the above RG equations to the t — t' model, we compute the flow starting with all the couplings 
set to the original value U of the on-site interaction. Under this initial condition, it is clear that Umtras = Uback and 
Uinters = Uumki all aloug the flow. The renormalized vertices show divergences at certain values of the frequency, 
measured in units of the cutoff Ec- We have 

= 1 + Uc/{7tH) \og{uj/E,) 

" l + Uc'/{nH) log(0./i?e) 

We interpret the divergences of the vertices in the same way as in the RPA, as signalling the development of an 
ordered phase in the system. 

The precise determination of the instability which dominates for given values of U and t' is accomplished by analyzing 
the respjCmse functions of the system. The procedure is similar to that followed in the study of one-dimensional electron 
systemsEJ. We adopt again a RG approach to their computation, which takes into account the scaling behavior of the 
interactions (Q) and (Q). We analyze ferromagnetic, antiferromagnetic, superconducting and CDW correlations. We 
assume that the actual instabilities of the model at low temperatures are some of these. 

It is easily seen that the operators related to charge-density-wave and s-wave superconducting instabilities do 
not develop divergent correlations at small oj. The phase diagram in the t' — U plane is drawn by looking at the 
competition among ferromagnetic, antiferromagnetic and d-wave superconducting instabilities. The ferromagnetic 
response function Rfm, for instance, is given by the correlation of the uniform magnetization, p^{Q^Lo) — ^^(O,^). 
The first perturbative terms for this object are built from a couple of one-loop particle-hole diagrams linked by 
the interaction. Each particle-hole bubble has a logarithmic dependence on the cutoff E^ with the prefactor c = 
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1/^1 - 4:{t'/ty. The iteration of bubbles can be taken into account by differentiating with respect to Ec and writing 
a self-consistent equation for Rfm, which turns out to be 

— 777^ = TTT^ i TT (Uintras + Uinters) -^RpM (9) 

OEc H^t Ec TT'^t Ec 

The antiferromagnetic response function Rafm can be dealt with in a similar fashion, by looking at correlations of 
the operator (QjCj) — p|(Q,cj). This leads to the RG equation 

ORafm _ 2c' 1 c' 1 

— T77^ ^ [^back + Uurnk) -f^^AFM (J-Uj 

OEc -K t Ec TT t Ec 

We recall that Uintras + Uinters and Uhack + Uumk havc the same flow, within the present model. Therefore, we may 
discern at once that whenever c > c' the ferromagnetic response function Rfm prevails over Rafm- 

Finally, there remains the response function for d-wave superconductivity Rscd, which is given by the correlation 
of the operator ^aj^|(k)aj^|(— k) — a^|(k)a^j^(— k) + h.c^ . The computation within the RG framework becomes 

now a little bit more subtle, since the diagrams at strictly zero total momentum display a log^ Ec dependence on the 
cutoff. As the derivative with respect to Ec is taken, it becomes clear that the logarithmic dependence left corresponds 
to the divergent density of states at the Van Hove singularity. The RG equation for Rscd reads then 

dRscd c \og{Ec/uj) _ logOBcM /^x 

dEc ~ 2nH Ec 2^2i^^™f- ^^rak) UsCd ^ij 

This equation also shows an homogeneous scaling of Rscd on uj/Ec, like in the previous cases. 

From inspection of Eq. (|l^), it is clear that divergent correlations in the d-wave channel arise for Uintras — Uumk < 0. 
According to the above results, this only happens for c < c', that is, outside the region of the phase diagram where 
Rfm > Rafm- This confirms that a ferromagnetic regime sets in for values of t' above the critical value t'c « 0.276<: 
at which c = c'. For values below t'c, there is a competition between Rafm and Rscd, which requires the analysis of 
the respective behaviors close to the critical frequency at which the response functions diverge. As a general trend, 
the response function Rscd dominates over Rafm in the regime of weak interaction, the strength being measured 
with regard to both the bare coupling constant and the value of the c' parameter. The reason for such behavior is that 
at weak interaction strength the RG flow has a longer run to reach the critical frequency, and at small frequencies 
the logarithmic density of states in Eq. makes Rscd to grow larger. The border where the crossover between 

the antiferromagnetic and the superconducting instability takes place is shown in the t' — U phase diagram of Fig. 
^. At sufhciently large values of U and small values of t' , the leading instability of the system turnS|-|ejiit, to be 
antiferromagnetism. This is in agreement with weak coupling RG analyses applied to the Hubbard modelBSy. 

Thus, there exists a region of the phase diagram where superconductivity is the leading instability. We remark 
that this result is obtained within a RG approach that provides a rigorous computational framework, with no other 
assumption than the weakness of the bare interaction. The instabilities are led by an unstable RG flow, and the 
singular behavior of the response functions is interpreted in the same fashion than in a standard RPA computation. 
The wide range for superconductivity is consistent with the results from quantum Monte Carlo corpoutationsEZl, as 
well as with results obtained by exact diagonalization of small clusters in the strong coupling regimet3. 

On the other hand, the renormalized interactions remain most part of the flow in the weak coupling regime for 
U /t < 1. We observe that the most interesting region in the phase diagram of Fig. ^ starts near U/t w 1, where the 
critical frequencies are ujc ~ 10~^Ec- This order of magnitude corresponds to sizeable critical temperatures 100 
K) if we assign to Ec a value of the order of the conduction bandwidth in the cuprates leV). 

The other relevant conclusion within our RG approach is the existence of a ferromagnetic regime in the t — t' 
model, above a certain value of the t' parameter. This is consistent with the results obtained in Ref. |l^ close to 
<' = 0.5t. Though our results refer to the weak coupling regime, they show that antiferromagnetic, ferromagnetic 
and superconducting phases are all realized in the t — t' Hubbard model. The superconducting instability has greater 
strength at the boundary with the antiferromagnetic instability, as it also happens in other approaches to high- Tc 
superconductivity^. 

In our case, the diagrams responsible for the apppearance of superconductivity cannot be interpreted in terms of the 
exchange of antiferromagnetic fluctuations, as in the work mentioned earlieiE3. Those diagrams which contain bubbles 
mediating an effective interaction between electron propagators are cancelled, to all orders, by vertex jserrections (see 
Fig. Superconductivity arises from the type of diagrams first studied by Kohn and Luttingeru. The strong 
anisotropy of tie Fermi surface greatly enhances the Kohn-Luttinger mechanism, with respect to its effect in an 
isotropic metalQ. 
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Our results support the idea that d-wave superconductivity and antiferromagnetism arise from the same type of 
interactions. Antiferromagnetism, however, does not favor the existence of superconductivity, but competes with 
it in the same region of parameter space. Similar physical processes seem to be responsible for the appearance of 
anisotropic superconductivity in systems of coupled repulsive ID chain£3. The Fermi surface of a single chain is 
unable to give rise to this type of superconductivity. A soon as this limitation is lifted, superconductivity occupies a 
large fraction of the phase diagram previously dominated by antiferromagnetic fluctuations. 

The results reported above can be extended to fillings away from the Van Hove singularity, provided that the 
distance of the chemical potential to the singularity is smaller than the energy scale at which the instability takes 
place. The chemical potential tends to be pinned to the singularity because of the nontrivial RG flow of the chemical 
potential itseliQ. Hence, these calculations can be applied to a flnite range of flUings around that appropriate to the 
singularity. 

In conclusion, we have shown that the t — t' Hubbard model at the Van Hove singularity exhibits a variety of 
instabilities at low energies or temperatures. The existence of these instabilities can be derived by RG methods which 
become exact at small couplings. We find that antiferromagnetism and d-wave superconductivity arise from the same 
interactions, and compete with each other in the same region of parameter space. 
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FIG. 1. Energy contour hnes about the Fermi level, with the Fermi hne passing by the saddle points A and B. 



FIG. 2. Different interaction terms arising from the flavor indexes A and B. 
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FIG. 3. Second order diagram renormalizing the different interactions in the model, with electron lines carrying flavor index 
A or B appropriate to each case. 

FIG. 4. Phase diagram in the {t', U) plane. The dotted lines are contour lines corresponding to the critical frequencies shown 
in the figure. 
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